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The effective interaction GXPFl for shell-model calculations in the full pf shell is tested in detail 
from various viewpoints such as binding energies, electro-magnetic moments and transitions, and 
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of shape coexistence suggest an advantage of Monte Carlo Shell Model over conventional calculations 
in cases where full-space calculations still remain too large to be practical. 
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I. INTRODUCTION 

The effective interaction is a key ingredient for the suc- 
cess of the nuclear shell model. Once a reliable interac- 
tion is obtained, we can describe various nuclear prop- 
erties accurately and systematically, which helps us to 
understand the underlying structure, and to make pre- 
dictions for unobserved properties. The pf shell for or- 
bitals lp3/2, lpi/2i 0/7/2 and O/5/2 is a region where the 
shell model can play an indispensable role, and is at the 
frontier of our computational abilities. In the pf shell 
one finds the interplay of collective and single-particle 
properties, both of which the shell model can describe 
within a unified framework. Since the protons and neu- 
trons occupy the same major shell, the proton-neutron 
interaction is relatively strong and one can study related 
collective effects such as T = pairing. The p/-shell 
nuclei are also of special interest from the viewpoint of 
astrophysics, such as the electron capture rate in super- 
novae explosions. For all these applications, a suitable 
effective interaction for p/-shell nuclei is required. 

Because of the spin-orbit splitting, there is a sizable 
energy gap between the /y /2 orbit and the other three or- 
bits (p3/2, Pi/2; fb/2)- Thus there exists an or Z = 28 
"magic" number inside the major shell with the oscillator 
quantum number Aqsc = 3. For shell-model calculations 
around this magic number, ^^Ni has often been assumed 
as an "inert" core. However, it has been shown that this 
core is rather soft0 and the closed-shell model for the 
magic number 28 provides a very limited description es- 
pecially for nuclei near N or Z=28 semi-magic. This 
"active-two-shell" problem is a challenge to both nuclear 
models and effective interactions. We shall discuss this 



point in this paper, and the word "cross-shell" refers to 
the N or Z=28 shell gap hereafter. Because of such cross- 
shell mixing, it is necessary to assume essentially the full 
set of pf configurations and the associated unified inter- 
action in order to describe the complete set of data and 
to have some predictive power. 

The effective interaction can in principle be derived 
from the free nucleon-nucleon interaction. In fact such 
microscopic interactions have been proposed for the pf 
shell H with certain success particularly in the be- 
ginning of the shell. These interactions, however, fail in 
cases of many valence nucleons, e.g., "^^Ca 0,3 and ^^Ni. 
Especially in the latter, the ground state is predicted to 
be significantly deformed in the full p/-shell calculation, 
contrary to its known double-magic structure. 

It has been shown that modifications in the monopole 
part of the microscopic interaction can greatly improves 
the description of experimental data. In fact, KB3 5] and 
KB3Gi interactions, which were obtained on the basis 
of the microscopic Kuo-Brown's G-matrix interaction^ 
with various monopole corrections, are remarkably suc- 
cessful for describing lighter p/-shell nuclei [A < 52). 
However, as we have pointed out^Tjj these interactions 
fail near ^^Ni. Therefore it is interesting to investigate 
to what extent the monopole modification is useful as a 
simple recipe for improvement of the microscopic inter- 
action and to what extent one has to go beyond this. 

One feasible way of modifying the microscopic inter- 
action for practical use is to carry out an empirical fit 
to a sufficiently large body of experimental energy data. 
In fact such a method has been successfully applied to 
lighter nuclei and has resulted in the "standard" eflfec- 
tive interactions of Cohen-Kurath |3| and USD for 
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the p and sd-shells, respectively. 

Along this line, we have recently developed a new effec- 
tive interaction called GXPFIJ],] for use in the pf shell. 
Starting from a microscopic interaction, a subset of the 
195 two-body matrix elements and four single-particle 
energies are determined by fitting to 699 energy data 
in the mass range A — A7 ^ 66. It has already been 
shown that GXPFl successfully describes energy levels 
of various p/-shell nuclei, such as the first 2+ states in 
even-even Z < 28 isotopes, low-lying states of 56,57,58]\jj^ 
and the systematics of the yrast l/2~, 3/2~ and 5/2~ 
states in Ni isotopes 0- Therefore it is now important 
to analyze the wave functions and examine some electro- 
magnetic properties to confirm further the reliability of 
the interaction. 

The aim of this paper is to present various results pre- 
dicted by GXPFl in order to clarify its applicability and 
limitation by comparing these results with experimental 
data. We also discuss what modifications of the micro- 
scopic interaction are needed in addition to the monopolc 
corrections in order to obtain a better description of nu- 
clear properties over the wide range of nuclei in the pf 
shell. 

In the derivation and previous tests of GXPFl, the 
Monte-Carlo Shell Model (MCSM) (lo| played a cru- 
cial role, since at that time it was the only feasible way to 
obtain shell-model eigenenergies for many states in the 
middle of the pf shell. On the other hand, in the present 
paper, most of the results have been obtained by the con- 
ventional Lanczos diagonalization method, which is now 
feasible with reasonable accuracy for most p/-shell nuclei 
owing to recent developments of an efficient shell-model 
code and fast computers. We will confirm our previous 
MCSM results by comparison with those by such conven- 
tional calculations, but there are still places where the 
MCSM is necessary. The pf shell is the current frontier 
of conventional methods, while the MCSM is applicable 
and useful in much larger model spaces, as demonstrated 
in Refs. [H mill El. 

This paper is organized as follows. In sectional the 
derivation of GXPFl is reviewed and its general proper- 
ties are analyzed in some detail. In section IIIII experi- 
mental data are compared with the results of large-scale 
shell-model calculations based on GXPFl and some its 
possible modifications. We analyze the structure of wave 
functions focusing on the core-excitations across the N 
or Z = 28 shell gap. A summary is given in section ITVI 



II. DERIVATION AND PROPERTIES OF THE 
EFFECTIVE INTERACTION 

In this section, we first sketch how we have derived the 
GXPFl interaction. The properties of GXPFl is then 
analyzed from various viewpoints. 



A. Derivation of GXPFl interaction 

An effective interaction for the pf shell can be specified 
uniquely in terms of interaction parameters consisting 
of four single-particle energies Ca and 195 two-body ma- 
trix elements V{abcd; JT), where a,b, ■ ■ ■ denote single- 
particle orbits, and JT stand for the spin-isospin quan- 
tum numbers. The eaS include kinetic energies. We take 
the traditional approach of evaluating the interaction en- 
ergy in zeroth-order perturbation theory for n nucleons 
outside of a closed shell for ^°Ca. We adjust the values of 
the interaction parameters so as to fit experimental bind- 
ing energies and energy levels. We outline the fitting pro- 
cedure here, while details can be found in 15]. For a set 
of N experimental energy we 
calculate the corresponding shell-model eigenvalues A^'s. 

We minimize the quantity — X]fcLi(^oxp ~ '^k)^ by 
varying the values of the interaction parameters. Since 
this minimization is a non-linear process with respect to 
the interaction parameters, we solve it in an iterative way 
with successive variations of those parameters followed by 
diagonalizations of the Hamiltonian until convergence. 

Experimental energies used for the fit are limited to 
those of ground and low-lying states. Therefore, cer- 
tain linear combinations (LC's) of interaction parameters 
are sensitive to those data and can be well determined, 
whereas the rest of the LC's are not. We then adopt 
the so-called LC method il6|, where the well-determined 
LC's are separated from the rest: starting from an initial 
interaction, well-determined LC's are optimized by the 
fit, while the other LC's are kept unchanged (fixed to the 
values given by the initial interaction). 

In order to obtain shell-model energies, both the con- 
ventional and MCSM calculations are used. Since we 
are dealing with global features of the low-lying spec- 
tra for essentially all p/-shell nuclei, we use a simplified 
version of MCSM: we search for a few (typically three) 
most important basis states (deformed Slater determi- 
nants) for each spin-parity, and diagonalize the Hamilto- 
nian matrix in the subspace spanned by these bases. The 
energy eigenvalues are improved by assuming an empir- 
ical correction formula, which is determined so that it 
reproduces the results of more accurate calculations for 
available cases. This method, the few-dimensional ap- 
proximation with empirical corrections (FDA*), actually 
yields a reasonable estimate of the energy eigenvalues 
with much shorter computer time. 

In the selection of experimental data for the fitting 
calculation, in order to eliminate intruder states from 
outside the present model space, we consider nuclei with 
A> 47 and Z < 32. As a result 699 data for binding and 
excitation energies (490 yrast, 198 yrare and 11 higher 
states) were taken from 87 nuclei : ■*^~^^Ca, ''^^^^Sc, 

47-52rpj 47-53, SSy 48-56q^ 50-58]y[j^ 52-60pg 54-61Qq 
56-66^1^ 58-63^^/60-642^; 62.64.65q; ^j^^ 64,65Gg_ 

assume an empirical mass dependence A~'^'^ of the two- 
body matrix elements similar to that used for the USD 
interaction 9], which is meant to take into account the 
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FIG. 1: Correlation of V{abcd; JT) between G and GXPFl. 
The matrix elements of T=0 and T=l are shown by open 
circles and crosses, respectively. For several matrix elements, 
corresponding quantum numbers are shown by using the no- 
tation 2ja2jt2jc2jd; JT. 



average mass dependence of a medium-range interac- 
tion jljj. We start from the microscopically derived 
effective interaction based on the Bonn-C potential 
which is simply denoted G hereafter. In the final fit, 70 
well-determined LC's are varied, and a new interaction, 
GXPFl, was obtained with an estimated rms error of 168 
keV within FDA*. The resultant single particle energies 
and two-body matrix elements are listed in Table 



B. Corrections to the microscopic interaction 

It is interesting to examine what changes have been 
made to the original G interaction by the empirical fit. 
Figure n shows a comparison between GXPFl and G for 
the 195 two-body matrix elements. One finds a strong 
correlation. On average, the T=0 {T=l) matrix elements 
are modified to be more attractive (repulsive). The most 
attractive matrix elements are T=0, with the largest ones 
belong to T=0 7575 in both GXPFl and G, as indicated 
in the figure, where the notation 7575 refers to the set of 
matrix elements V(abab; JT) with a = /7/2 and 6 = /5/2- 
It was stressed in [1^13 that this strong T—0 interaction 
between the orbits j> = 1 + 1/2 and j< = Z — 1/2 is 
an important feature in all mass regions including the 
sd and p shells. In the p shell variation of these type of 
matrix elements relative to the Cohen-Kurath interaction 
leads to improvement 0,1131 • One can see in Fig. 
that the the G and fitted (GXPFl) values for the j>-j< 
(7575) interactions are very similar. 



There are seven matrix elements in which the difference 
between GXPFl and G is greater than 500 keV. These 
matrix elements are listed in Table Ull It is remarkable 
that these largely-modified matrix elements are either of 
the diagonal V{abab; JT) type which contributes to the 
monopole corrections or of the monopole pairing (J — 0, 
T ~ 1) type. As for the former type, the 7373 and 7575, 
T — matrix elements with large J are modified to be 
more attractive, contrasting to the relatively small cor- 
rections for small J matrix elements. The two monopole 
pairing matrix elements shown in the table are both re- 
lated to the /5/2 orbit. As a result of the empirical fit, 
the pairing between P3/2 and /5/2 (3355) is modified to 
be strongly attractive, while that between and /5/2 
(7755) is made to be less attractive. 



C. Monopole properties 

In order to investigate basic properties of an effective 
Hamiltonian, it is convenient to decompose it into the 
monopole part and the multipole part 21] as. 



H — Hn 



H 



M- 



(1) 



The monopole part Hm plays a key role for describing 
bulk properties such as binding energies and shell-gaps 
[1| , since it determines the average energy of eigenstates 
in a given configuration. The monopole Hamiltonian is 
specified by the angular-momentum averaged two-body 
matrix elements: 



V{ab;T) 



J2j{'2J+l)V{abab; JT) 
Ej(2J+1) ' 



(2) 



where the summations run over all Pauli-allowed values 
of the angular momentum J. 

Figure|21shows the matrix elements V{ab; T). As a ref- 
erence, we consider also the KB3G P| interaction. Since 
this interaction gives an excellent description for A < 52 
nuclei, it is expected that the monopole matrix elements 
of GXPFl are similar to those of KB3G at least for those 
involving f7/2- In fact, both T — and T = 1 average 
matrix elements of the f7f7, f7p3, f7f5 and f7pl orbit pairs 
are rather close between GXPFl and KB3G. The T =1 
matrix elements for p3p3, p3f5 and f5pl are also not very 
different, which are important for describing Z < 28, 

> 28 nuclei. On the other hand, the similarity is lost 
in other matrix elements, especially for T — matrix el- 
ements between ^3/2,1/2 orbits. Therefore, GXPFl and 
KB3G could give a very different description for nuclei 
with Z,N > 28. 

In order to confirm this observation, we have carried 
out the FDA* calculations by using KB3G for the similar 
set of nuclei included in the fitting calculations. Tabic UTTI 
summarize the estimated rms difference between calcu- 
lated energies and the experimental data for both GXPFl 
and KB3G. The data are classified into four groups from 
the viewpoint of the location in the isotope table : (a) 
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TABLE I: Two-body matrix elements V{abcd; JT) (in MeV) of GXPFl interaction. Single particle energies are taken to be 
—8.6240, —5.6793, —1.3829 and —4.1370 MeV for the /7/2, P3/2, /5/2 and pi/2 orbit, respectively. For calculations of mass A 
nuclei, these two-body matrix elements should be multiplied by a factor . 
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Z,N < 28, (b) Z or TV = 28, (c) Z < 28, N > 28 and (d) 
Z,N > 28. The rms deviations are estimated for each 
group, and the yrast and yrare states are considered sep- 
arately. 

It is clearly seen that the rms deviations by GXPFl 
are almost the same in all groups. On the other hand. 



those by KB3G become larger for the group (b) and (d) 
in comparison to (a) and (c). The deviation is especially 
large for the group (b), where the major differences can 
be foimd in "''^Ni, '^^Co and '"'^Ni. Since the property 
of the core-excitations appear directly in the low-lying 
spectra in these semi-magic nuclei, this result suggests 
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TABLE 11: Comparison of the two-body matrix elements 
V{abcd; JT) (MeV) {A = 42) for which the difference between 
G and GXPFl is larger than 500 keV. 



2ja 


2jb 


2jc 


2id 


J 


T 


G 


GXPFl 


difference 
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7 
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-2.2033 
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-0.7637 
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+0.6043 



0.5 



TABLE IIL Comparison of rms deviations (MeV) between the 
experimental excitation energies and those calculated from 
GXPFl and KB3G. Theoretical energies are estimated by the 
FDA*. The numbers in the parentheses show the number of 
data included in the calculation. 



group 


state 


GXPFl 


KB3G 


(a) ^, < 28 


yrast 


0.154(136) 


0.235(129) 




yrare 


0.201(45) 


0.263(23) 


(b) Z OY N ^28 


yrast 


0.184(92) 


0.647(87) 




yrare 


0.195(57) 


0.802(44) 


(c) ^ < 28, A > 28 


yrast 


0.145(129) 


0.296(126) 




yrare 


0.145(75) 


0.302(55) 


(d) ^, A" > 28 


yrast 


0.186(55) 


0.401(51) 




yrare 


0.187(23) 


0.458(23) 




f7f7 f7p3 f7f5 f7p1 p3p3 p3f5 p3p1 f5f5 f5p1 p1p1 



FIG. 2: Comparison of the monopole matrix elements 
V{ab- T) {A = 42) between G, GXPFl and KB3G, which are 
shown by circles, squares and diamonds, respectively. Lines 
are drawn to guide the eyes. The orbit-pair label "f5p3" 
stands for a = f-j/2 and b = P3/2, for example. 



that the core-excitations are not well described by the 
KB3G interaction. 

The rms deviations are in general larger for yrare states 
than the yrast states, especially in KB3G. This also sug- 
gests the failure in the description of the core-excitation, 
which is expected to appear more directly in the yrare 
states than the yrast states. Another possible reason is 
that the FDA* becomes less accurate for yrare states. 

From this table we cannot infer that the description 
by GXPFl is better than KB3G for group (a) and (c), 
because the rms deviations shown in the table are not ex- 
act values but the results of the FDA* with typical accu- 
racy of about 200 keV. Since GXPFl is determined from 
FDA*, the rms deviations of GXPFl is naturally smaller 
than those by other interactions within the FDA*. In 
fact, in light p/-shell nuclei with more precise computa- 
tions, we can find several examples where KB3G gives a 
better description than GXPFl. 



D. Collective properties 

It has been pointed out that the multipole part 
of the Hamiltonian Hm is dominated by several terms, 
such as the pairing and the quadrupole-quadrupole in- 
teractions, which determine collective properties of the 
effective interaction. It is useful to investigate these col- 
lective aspects of GXPFl and compare with that of other 
interactions. 

According to the prescription in Ref. |23|, can be 
expressed in both the particle-particle (p-p) representa- 
tion and the particle-hole (p-h) representation. In the 
p-p representation, Hm is expressed in terms of particle- 
pair bases [cjjcj]'^"^, in which two nucleons are coupled to 
the good angular momentum J and the isospin T. The 
coefficient matrix is diagonalized for each JT channel. 
By using the resultant eigenvalues iS^"^, we obtain 

Hm = Ei^PjTaPjTa + (onc-body terms), (3) 

JTa 

where Pjj'a denote particle-pair creation operators which 
are linear combinations of the particle-pair bases, and 
their structure is determined by the corresponding eigen- 
vectors. (Here, for simplicity, we omit the summations 
over z-components of both spin and isospin.) The most 
important contributions come from the JT = 10, 20 and 
01 terms with large negative eigenvalues E'^'^ , which cor- 
respond to the T = 0, small-J pairing and the usual 
T = 1 monopole pairing, respectively. 

Similarly, in the p-h representation, by using the 
particle- hole (density) bases [cjjCf,]'^'^ with spin- isospin 
quantum numbers At, the diagonalization of the coeffi- 
cient matrix determines the structure of multipole oper- 
ators Qxrk and the corresponding strengths e^'^, leading 
to an expression 

Hm ~ 'Y^^k^QxTk ■ Qxrk- (4) 

Arfc 

The symbol • stands for a scalar product with respect 
to both spin and isospin. (We adopt here a slightly 
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TABLE IV: Comparison of the collective strengths (MeV) 
between G, GXPFl and KB3G. The mass number A = 42 is 
assumed. 



Interaction 


^01 




^20 






e" 


G 


-4.20 


-5.61 


-2.96 


-3.33 


-1.30 


+2.70 


GXPFl 


-4.18 


-5.07 


-2.85 


-2.92 


-1.39 


+2.67 


KB3G 


-4.75 


-4.46 


-2.55 


-2.79 


-1.39 


+2.47 



different definition of e^"^ from that in Ref. |^ by 
a pliase factor (—1)^+'^.) Large eigenvalues appear for 
At = 20, 40 and 11, which are interpreted as the usual 
isoscalar quadrupole, hexadecapole and Gamow- Teller 
interactions, respectively. 

In Table IIVI the collective strengths, i.e. the largest 
(smallest) eigenvalues E'^'^ or e^'^ in each spin-isospin 
channel, are shown for G, GXPFl and KB3G. According 
to the comparison between G and GXPFl, in general, the 
empirical fit has reduced these strengths by at most 12% 
except for the At = 40 term. Especially, the reductions 
in the T = terms E^'^ and are relatively large, 
contrary to the general observation that T — matrix 
elements are on average modified to be more attractive 
by the empirical fit (see Fig. QJ. This means that the 
attractive modification has been applied mainly to the 
monopole terms (which do not enter into the numbers of 
Table eg. 

For the comparison between GXPFl and KB3G, it 
can be seen that the T = (T = 1) pairing strength 
is stronger (weaker) in GXPFl. Such a difference is 
expected to affect the description of the structure for 
Z N nuclei. In GXPFl, E^° is larger than E°^ by 
about 20%, which is consistent with an estimate by the 
mean-field calculations using the standard seniority 
pairing. It has also been shown that the T — 
pairing strength is larger than T = 1 in the density- 
dependent Gogny force. Note that the situation is oppo- 
site in KB3G. On the other hand, the multipole strengths 
e^^ are very similar for GXPFl and KB3G, although the 
strengths of GXPFl are slightly larger. 

By using e^"^ , we can evaluate the collective 
quadrupole-quadrupole (QQ) strength between like- 
nucleons (proton-proton or neutron-neutron) and that 
between protons and neutrons. For GXPFl, the strength 
of Qp ■ Qp or Qn • Qn is —0.96 MeV, while that of Qp ■ Qn 
is —7.82 MeV. This result shows the dominance of the 
proton-neutron part in the collective QQ interaction, as 
in heavier nuclei 0, IH, 0, 113 • Such Qp-Qn dominance 
can be seen in other interactions, although the ratio of 
Qp-Qn to Qp-Qp (or Qn-Qn) is different (8.1 for GXPFl, 
6.8 for G, 7.1 for KB3G). 



E. Spin-tensor decomposition 

In order to analyze the structure of an effective interac- 
tion, the spin-tensor decomposition llTj is useful, since it 



gives physical insights from a different viewpoint. In the 
following discussions, we consider the two-body interac- 
tion Vm in the multipole part Hm- Therefore the results 
are free from the monopole effects. We again take KB3G 
as a reference interaction. It is essentially the same as the 
Kuo-Brown's renormalized G-matrix interaction [2| after 
the monopole subtraction. (All of the results in this sec- 
tion are for the matrix elements evaluated at A = 42.) 
An overview of the correlation in the Vm between these 
interactions are shown in Fig. |21 It can be seen that 
the correction to the microscopic interaction imposed by 
an empirical fit contains sizable non-monopole compo- 
nents. In addition, the microscopic interactions, G and 
KB3G (^ KB) are not identical even in the multipole 
parts, as seen in the lower part of Fig|3| The difference 
between G and KB3G looks not necessarily smaller than 
that between G and GXPFl, which also suggests that 
the present correction to G is in a reasonable magnitude. 

We first transform the jj-coupled two-body matrix el- 
ements {jajbJT\VM\jcjdJT) into the LS'-coupled form 
{lalbLSJT\VM\lJdL'S'JT), then carry out the spin- 
tensor decomposition of the two-body interaction as 

VM = Y.Vk=Y.U'^-X\ (5) 

k k 

where the operators U*' and X*' are irreducible tensors 
of rank k in the space and spin coordinates, respectively. 
The interaction components Vk represent the central {k = 
0), spin-orbit {k = 1) and tensor (fc — 2) parts. The spin- 
orbit part includes both the normal part [S = S' = \) 
and the anti-symmetric spin-orbit {S ^ S') part. 

Figure^lshows the central components. For the T — 
matrix elements, G and GXPFl are very similar to each 
other. Sizable modifications of G by the empirical fit can 
only be seen in the {lalblJd L)=(3333 0) and (3131 2) 
components, which are made more repulsive by about 
0.4 and 0.3 MeV, respectively. Other matrix elements 
of GXPFl are very close to those of G. On the other 
hand, GXPFl deviates from G in various T = 1 matrix 
elements. The modifications to G are in the repulsive 
direction for the matrix elements (3333 0) and (3333 2) 
which are related to monopole and quadrupole pairing 
in the / orbit, respectively, while the (3311 0) matrix 
element is made to be more attractive, which corresponds 
to the monopole pairing between / and p orbits. 

For KB3G, it is remarkable that T = 0, S" = 1 matrix 
elements are very close to those of GXPFl, including the 
(3333 0) matrix element which deviates from G signifi- 
cantly. Considering that the origin of KB3G and GXPFl 
is very different, this similarity is surprising. However, 
there are large differences in the T — 0, S = matrix 
elements, where the absolute values of the KB3G matrix 
elements are smaller in most cases than those of GXPFl 
(and also G). The only exception is the (3333 1) matrix 
element. On the other hand, most of the T — 1 matrix 
elements of KB3G are very close to those of G rather 
than GXPFl for both S = and 1. One can only find 
small deviations from G in the (3333 0) and (3333 6) 
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FIG. 3: Correlation of monopole-free two-body matrix ele- 
ments Vai between (upper) GXPFl and G, (middle) GXPFl 
and KB3G and (lower) G and KB3G. 
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FIG. 4: Comparison of the central T = (lower panel) and 
T — 1 (upper panel) components. Results for GXPFl, KB3G 
and G are shown by solid lines, dotted lines and circles, re- 
spectively. The quantum numbers lahlcldLL' SS' of the LS- 
coupled matrix elements {lalbLSJT\V\lcldL' S' JT) are shown 
along the horizontal axis. 



elements, where the former is more attractive and the 
latter is more repulsive. This similarity between G and 
KB3G indicates that the T = 1 central part is converged 
in these two different G-matrix calculations. 

In Fig. [S] the tensor components are compared. The 
T — matrix elements are relatively large for AL = 
2. The corrections to G are relatively large for the 
{UbUd LL')=(3331 02), (3331 42), (3311 02), (3131 42) 
and (3111 42) matrix elements, which are all in the repul- 
sive direction, while those in the attractive direction are 
found in (3333 64) and (3331 64) matrix elements. These 
corrections are at most 0.2 MeV. The T = 1 components 
are in general very small compared to T = 0. This is due 
to the fact that T = tensor interaction is dominated 
by a matrix element between two nucleons with Z = to 
I = 2 whereas T = 1 must have odd /. 

The G and KB3G values show rather large differences 
in cases with largest magnitudes, while the GXPFl fa- 
vors the G values with small corrections. Namely, in the 
T = cases, (3333 20), (3331 02) and (3111 20) matrix 
elements, the attraction of KB3G is much weaker, rather 
outstandingly, than that of G (and also GXPFl). 

The normal spin-orbit components are shown in Fig. 
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FIG. 5: Comparison of the tensor components, 
are the same as in Fig. 0] 



Conventions FIG. 6: Comparison of the spin-orbit components. Conven- 
tions are the same as in Fig. 0| 



For the T = elements, relatively large attractive cor- 
rections to G are found in the {lahlJ-d ii')=(3333 22), 
(3333 44) and (3131 44) matrix elements. Such correc- 
tions to the latter two matrix elements are not found in 
KB3G. Similarly, large attractive corrections exist also in 
the T = 1, (3331 12), (3331 32), (3331 54) and (3131 33) 
matrix elements, all of which are absent in KB3G. It can 
be seen that the spin-orbit matrix elements of KB3G and 
G are very close in both cases of T = and 1. Thus, we 
infer that there may be contributions to the spin-orbit 
interaction that are not present in the two-nucleon G 
matrix, perhaps from an effective three-nucleon Hamil- 
tonian. 

Figure shows the anti-symmetric spin-orbit compo- 
nents. In general, matrix elements of G are small as in the 
case of the normal spin-orbit components. Nevertheless, 
one can again see nearly perfect similarity between G and 
KB3G in almost all matrix elements. The empirical fit 
has resulted in several large corrections to G which are in- 
consistent with KB3G, such as T = 0, {lahlJd LL' SS')~ 
(3331 32 01), (3331 54 01), (3331 44 10) and T = 1 
(3333 10 10), (3311 10 10), (3111 32 10). Such relatively 
large matrix elements appear also in_other empirical ef- 
fective interactions such as FPMI3 261 and TBLC8 H^. 
Further investigation is needed to find if these deviations 
are significant in terms of the errors inherent in the fits 



to data. 

Summarizing the results of the spin-tensor decompo- 
sition of Vm , there are overall, reasonably good similar- 
ities between G and GXPFl in relatively large matrix 
elements such as the central and the T = tensor com- 
ponents. The corrections due to the empirical fit become 
sizable in several specific matrix elements, especially for 
T = 1. On the other hand, although many of matrix 
elements of KB3G are very close to those of G, one sees 
that several T = 5 = central and T = tensor matrix 
elements are rather different. This suggests significant 
changes in microscopic calculations of the effective inter- 
action or their input from the nucleon-nucleon interac- 
tion that have evolved from the original calculations of 
Kuo-Brown to the more recent results of Ref . l3| . 



F. Monopole fit 

According to previous subsections, the empirical fit 
gives rise to sizable corrections to the microscopic inter- 
actions mostly in the monopole part and several specific 
matrix elements such as the monopole pairing. There- 
fore we come to a natural question: to what extent can 
we improve the microscopic interaction for practical use 
with more restrictive corrections. In order to assess this 
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FIG. 7: Comparison of the anti-symmetric spin-orbit compo- 
nents. Conventions are the same as in Fig. 2] 



approach, we have carried out another fit by varying only 
the monopole parts, and the monopole-pairing {T — I, 
J — 0) and quadrupole-pairing (T = 1, J — 2) ma- 
trix elements. The number of parameters are 20, 10 and 
36 for these components, respectively. Thus in total 70 
parameters are varied including 4 single particle ener- 
gies. Note that the T = components are varied only 
in the monopole part. The mass dependence A~^'^ is 
also assumed. In the final fit, 45 best-determined LC's 
are taken, and the resultant interaction is referred to as 
GXPFM. The estimated rms deviation is 226 keV within 
FDA* for 623 energy data. Although this number sounds 
quite successful, GXPFM fails to describe N or Z—28 
semi-magic nuclei. In fact, for these nuclei, the estimated 
rms deviation increases to 267 keV for yrast states (87 
data) and to 324 keV for yrare states (44 data). 

We first consider the multipole part of the Hamilto- 
nian. The collective strengths of the p-h interactions of 
GXPFM are found to have reasonable values. For exam- 
ple, e^o = -2.85, = -1.44 and e" = -^2.59 (MeV), 
which are quite similar to those of GXPFl (see Table 
II V|) . As for the p-p channel, since the multipole part of 
GXPFM is different from G only for the T = 1 compo- 
nents, the strengths and are the same as those of 
G. On the other hand, the monopole pairing strength is 
reduced to E'^-^ = —3.88, which is much weaker than that 
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FIG. 8: Comparison of the monopole matrix elements 
V(ah; T) for the GXPFl and GXPFM interactions, which are 
shown by squares and triangles, respectively. Conventions are 
the same as in Fig. |21 



of GXPFl and will compensate the large T = pairing 
strength. 

Figure |S1 compares the monopole matrix elements 
y(a6; T) for GXPFl and GXPFM. The /7/2-related ma- 
trix elements of GXPFM take similar values to those of 
GXPFl (and other interactions), although the T = 0, 
f7f5 matrix element of GXPFM is more attractive than 
that of GXPFl by 300 keV. Note that the monopole 
Hamiltonian is specified by the linear combinations 
\V(ah; 1) -f \V{ah; 0) and V(ah; 1) - V(ah; 0), where the 
former determine, roughly speaking, the mass depen- 
dence and the latter affects the isospin dependence. The 
first combination is extremely well determined and al- 
most identical between GXPFl and GXPFM (differences 
are less than 50 keV). 

Relatively large differences can be seen in the T = 
1, p3p3, f5f5 and plpl matrix elements, which are 
all related to the monopole pairing. The former two 
of GXPFM are less attractive than those of GXPFl, 
which is required to keep the collective monopole pair- 
ing strength weaker than that of GXPFl, as mentioned 
above. For the T = matrix elements, the p3f5 and 
p3pl components are very different between GXPFl and 
GXPFM. Although the structure of Ga and Ge isotopes 
is rather sensitive to these matrix elements, most of those 
isotopes cannot be included into the fit because iV is too 
close to 40. This point may be clarified better in a future 
study including higher orbits. 



III. RESULTS AND DISCUSSIONS 

In this section, the results obtained from large-scale 
shell-model calculations with GXPFl are given and com- 
pared with available experimental data. In our previous 
studies of p/-shell nuclei 0, IH, H^, we took ad- 
vantage of the MCSM, which was the only feasible way 
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to evaluate the eigenvalues of the shell-model Hamilto- 
nian in the middle of the pf shell. In the present study, 
most of shell-model calculations are carried out in a con- 
ventional way by using the shell-model code MSHELL 
[sof. MSHELL enables calculations with M-scheme di- 
mensions of up to ^ 10^. With this capability we can 
handle essentially all low-lying states in the p/-shell nu- 
clei with a minimal truncation of the model space. We 
can also confirm the reliability of the previous MCSM cal- 
culations by comparing the both results. On the other 
hand, for non-yrast states in the middle of the pf shell, 
the MCSM it is still necessary. 

In the following discussions, the truncation order t de- 
notes the maximum number of nucleons which are al- 
lowed to be excited from the /7/2 orbit to higher three 
orbits P3/2, /5/2 and Pi/2, relative to the lowest filling 
configuration. The latter three orbits are expressed sim- 
ply as r hereafter. Most results for Ca, Sc, Ti, V, Cr are 
exact (no truncation) while t > 6 for Mn, and t > 5 or 
better for other isotopes. 

In this paper we fully cover the data on magnetic and 
quadrupole moments. Electro-magnetic transitions are 
discussed for some representative nuclei. Gamow- Teller 
beta decay will be covered in subsequent work. Electro- 
magnetic transition matrix elements are calculated by 
using the effective g-factors and the effective charges 
adopted in sections IlII CI and IIII Dl respectively. 



A. Closed core properties 

First we consider the role of the Z or iV = 28 closed 
shell. This is important for the unified shell-model de- 
scription of p/-shell nuclei, since it is convenient to base 
a truncation scheme on the stability of such a closed 
core. It has been shown jlj] that the ^^Ni core is rather 
soft in comparison to the '^^Ca core. More precisely, the 
probability of the (77/2)^^ configuration in the ground- 
state wave function of ^^Ni is much smaller than that of 
the (/7/2)^ configuration in ''^Ca. It is attributed to the 
strong proton-neutron interaction. The evolution of such 
a closed core in various isotope chains is of interest. 

Hereafter, a group of configurations in which both pro- 
ton and neutron orbits are maximally filled will be 
denoted by the ^^Ni closed-shell configuration. For each 
isotope with Zy valence protons and valence neu- 
trons on top of the ^"Ca core, such configurations are 
those given by 7r(/7/2)™(r)^"-™i/(/7/2)"(r)^''-" with 
m = min{Z„, 8} and n = min {A^„,8}. FigureElshows the 
probability of the closed-shell configurations in the cal- 
culated ground-state wave function. The lowest T = 
states are considered for odd-odd N = Z isotopes, al- 
though they are not necessarily the ground states. The 
truncation order is taken to be sufficiently large for de- 
scribing these states. 

The probability takes the minimum value (~40%) 
around *^Cr, which can be understood as a result of 
the large deformation. In the Ti isotopes, the probabil- 
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FIG. 9: The probability of ^^Ni closed-shell configurations in 
the calculated ground-state wave functions. 

ity becomes smallest around iV 24 and then increases 
monotonously for larger TV. On the other hand, in the 
case of Fe and Ni, one can see a local maximum (~60%) 
at iV = 28 and the second minimum (~45%) around 
TV = 32 ~ 34. Thus it is not justified to assume an in- 
ert ^^Ni core even for Ni isotopes. As for Zn isotopes, 
the probability shows rather smooth behavior with mod- 
erate values ~65%, suggesting that the effects of core- 
excitations are similar over the wide range of N . 



B. Binding energy 

Binding energies are obtained by adding suitable 
Coulomb energies to the shell-model total energies. In 
the present study the Coulomb energies are evaluated by 
using an empirical formula: 

= Ktt '^*''^^ - V-KvT^V + e^TT, (6) 

where tt and v denote the number of valence protons 
and neutrons, respectively. The adopted parameters are 
14^ = 0.264, V;^ = 0.038, and 7.458 (MeV). The 
same form was adopted in Ref. for describing light 
p/-shell nuclei with a different parameter set. In the 
present approach, since we consider a wider mass region, 
the parameters are determined by fitting to the energy 
difference between 36 pairs of isobaric analog states with 
masses A = 47 ~ 74. In Ref. '35|, the binding energies 
are also studied systematically by using the KB3 interac- 
tion for many /7/2 shell nuclei up to ^^Ni, and quite sim- 
ilar values of these parameters are proposed. They have 
attained an rms deviation of 215 keV between theory and 
experiment for 70 nuclei of ^ = 42 ~ 56. However, the 
discrepancy in ^^Ni (overbinding) is significantly large in 
comparison to that of neighboring nuclei. 

In Fig. ^1 the deviation of calculated binding ener- 
gies from experimental values are shown for each isotope 
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chain as a function of the neutron number N. The the- 
oretical values are obtained by GXPFl. For the trun- 
cations, all results are exact for Ca, Sc, Ti, V and Cr 
isotopes. For other isotopes, the results are obtained in 
a subspace of t = 5 or larger. Judging from the con- 
vergence pattern as a function of t, the underbinding 
due to this truncation is typically smaller than 200 keV. 
Note that the mass range included in the figure is much 
wider than that used in the fitting calculations for de- 
riving GXPFl, and covers regions of nuclei which may 
not be well described just by the pf shell. It can be 
seen that the agreement between theory and experiment 
is quite good over the whole mass range included in the 
fit {A — A7 — 65) and is also reasonable for many nuclei 
which were not included in the fit. 

In the neutron-rich side, we can find relatively large 
deviations. In Ca and Sc isotopes, the calculations give 
overbinding at the end of the isotope chain, iV = 32 and 
34, respectively, with large error-bars in the experimental 
data. Note that, for ^^Ca, the experimental data was 
taken from Ref. [S^l, where the value Qa_=7.9(5) MeV 
is adopted. On the other hand, in Ref. |34l| . the measured 
value Qp-=5.7{2) MeV is presented. The prediction by 
GXPFl is 5.9 MeV, in good agreement with the latter 
value. It is important to have improved experimental 
data for this neutron-rich region of nuclei. For the Ti 
isotopes, the agreement is satisfactory along the whole 
isotope chain. In other isotopes, the deviation between 
theory and experiment becomes sizable around A'^ 35, 
and, for larger N, the calculation shows a systematic 
underbinding. For a fixed N, the deviation is largest for 
Cr, and it is moderate for V, Mn, Fe, Ge, and it is small 
for Co, Ni, Cu, Zn, Ga. We relate this to an increasing 
importance of the neutron gg^2 orbit which should have 
its maximal effect in the middle of proton /7/2 shell (Cr) 
due to deformation. 



C. Magnetic dipole moments 

The magnetic dipole moments are calculated and com- 
pared with experimental data in Table FVl The magnetic 
moment operator used in the present calculation is 



9il, 



(7) 



where gs and gi are the spin and the orbital g-factors, 
respectively. By using the free g-factors gs = 5.586, gi — 
1 for protons and gs — —3.826, gi = for neutrons, 
the agreement between the calculation (/iti°°) ^nd the 
experiment (/ioxp) is in general quite good. 

Figure lTTl shows the comparison of magnetic dipole mo- 
ments between the experimental data and the shell-model 
predictions. The agreement appears to be good except 
for several cases to be discussed. As in the case of the 
sd-shell 0, the description is successful already by us- 
ing the free-nucleon g-factors. It is well known that ex- 
perimental magnetic moments deviate strongly from the 
single-particle (Schmidt) values, as seen for example on 
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FIG. 10: Comparison of experimental binding energies with 
the shell-model energies. The upper (lower) panel shows the 
results of even- (odd-) Z isotope chains. Data are taken from 
Ref. lal. 



the left-hand side Fig. 37 of |33 • For the sd and pf shells 
we find that the configuration mixing within the shells is 
enough to fully account for observed magnetic moments. 
This is in contrast to Gamow- Teller beta decay where the 
matrix elements of the isovector spin factor, gs, are sys- 
tematically reduced by factors of 0.77 in the sd shell 
and 0.74 for the pf shell 36] . Thus, the good agreement 
obtained for the magnetic moments with the free-nucleon 
operator is interpreted as cancellation of the quenching 
observed in Gamow- Teller decay with enhancements in 
the spin and orbital electromagnetic operators due to ex- 
change currents 
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TABLE V: Comparison of experimental magnetic dipole mo- 
ments ficxp (in unit ^n) with theoretical values f/^-^° and fi^^, 
which are calculated by using the free and the effective g- 
factors, respectively. Most of the data are taken from Ref. 

except for from Ref. [H, ^ from Ref. [s^, " from 
Ref.*!^, from Ref. " from Ref. li^ and ^ from Ref. 

|43J." All results for Ca, Sc, Ti, V, Cr are exact, while Z; > 6 
for Mn, and t > 5 for other isotopes. 
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0, 


.492 




3/2- 


0, 


.185 


+0.50(6) 


0, 


.679 


0, 


.543 


«8Zn 


2+ 


1, 


.077 


+0.87(9)' 


1, 


.489 


1, 


.363 


™Zn 


2+ 


0, 


.885 


+0.76(8)' 


3, 


.757 


3, 


.704 


«6Ga 


(2) + 


0, 


.066 


1.011(18) 


1, 


.056 


0, 


.874 




(9+) 


3, 


.043 


4.2(9) 


4, 


.443 


4, 


.422 


«^Ga 


3/2- 


0, 


.000 


+1.8507(3) 


1, 


.792 


1, 


.786 


6«Ga 


5/2- 


0, 


.359 


1.40(65) 


1, 


.310 


1, 


.584 


1+ 


0, 


.000 


0.01175(6) 


0, 


.188 


0, 


.003 


'^''Ga 


3/2- 


0, 


.000 


+2.01659(4) 


1, 


.850 


1, 


.840 


^^Ga 


3/2- 


0, 


.000 


+2.56227(2) 


2, 


.782 


2, 


.748 


69 /"I 

Ge 


5/2" 


0, 


.000 


0.735(7) 


1, 


.063 


0, 


.757 


™Ge 


2+ 


1, 


.039 


+0.936(52) 


0, 


.789 


0, 


.671 


^^Ge 


1/2- 


0, 


.000 


+0.547(5) 


0, 


.346 


0, 


.250 


^2Ge 


5/2- 


0, 


.175 


+1.018(10) 


1, 


.255 


0, 


.909 


2+ 


0, 


.834 


+0.798(66) 


2, 


.314 


2, 


.469 



However, for the Ni and Zn isotopes, calculated values 
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FIG. 11: Comparison of experimental magnetic dipole mo- 
ments with the sheU-model results, which were obtained by 
using the free (upper panel) and effective (lower panel) nu- 
cleon g- factors. All data in Table ^ are included for which 
the sign is measured experimentally. 



are systematically larger than the experimental ones by 
typically 0.2 ^ O.S/iat. Such difficulties can be remedied 
to a certain extent by introducing the effective g-factors, 
as shown in the same table (/it^)- present cal- 

culation, we took g";^ = 0.9(7g''°°, gi = 1.1 and —0.1 for 
protons and neutrons, respectively, which were chosen 
from an estimate by the least-squares fit. It can be seen 
that the description of odd-A nuclei is systematically im- 
proved. 

Nevertheless, deviations on order of 0.3/iAr remain for 
the 2+ state of Zn nuclei. We have examined the effect 
of a more general effective Ml operator which contains 
the [y2'S]'-^-' term, but it turns out that the deviations 
for the Zn nuclei could not be remedied. It may be re- 
quired to include the 59/2 orbit to improve the descrip- 
tion of these states, as also discussed in Ref. One 



may also want to reexamine the systematic uncertainties 
which may exist in the effective transient fields which are 
used to deduce magnetic moments from the experimental 
data. 

There are a few other cases where we find large differ- 
ences between theory and experiment. Several of these 
can be interpreted as a consequence of incorrect mixing 
of two closely-lying states. For example, in '"^^Fe, the first 
and the second 7/2~ states are separated only by 92 keV 
experimentally. The GXPFl predicts these states in the 
reversed order. It is also the case for two 3/2~ states 
of ^^Zn, where the experimental energy separation is 92 
keV. 

We also find notable differences between theory and 
experiment for ^OCr 8+, ^^Mn 5+, ^^Co 3+, ''^Ti 7/2" 
and '^"'^Cr 3/2^ where experimental uncertainties are also 
large. More precise measurements of these are required. 
The deviations in ^^Cu 5/2~ and ^^Cu 1+ may be at- 
tributed to the effect of ^9/2 orbit. The deviations in 
^"Zn and ^^Ge are naturally understood as a result of 
the insufficient model space. 

D. Electric quadrupole moments 

The electric quadrupole moments are given in Table 
IVII The eff'ective charges e^r = 1.5, = 0.5 are adopted 
in the present calculations. The correlation between the 
calculated electric quadrupole moments and the experi- 
mental data is shown in Fig. 1121 One can find that, in 
general, the deviation of the theoretical prediction from 
the experimental data becomes large where the exper- 
imental error bar is also large, except for a few cases. 
The sign of the calculated Q{3/2~) for ^^Co is opposite 
to the experimental data. This result suggests an incor- 
rect mixing of the first and the second 3/2^ states in 
the calculation, which are close in energy (380 keV ex- 
perimentally). The calculated 3/2^ state lies 589 keV 
above 3/2]^ state with Q^-f 0.138eb, which is consistent 
with this interpretation. The description of ^*Zn 2+ is 
also very poor, but the uncertainty in the data is large. 
The description is unsuccessful for ''"Zn, ^Ge and ""^Ge, 
indicating the need for introducing the 39/2 orbit. 

E. Systematics of 2^ states 

The first 2+ state of an even-even nucleus is a good sys- 
tematic measure of the structure. The left panel of Fig. 
1131 shows the excitation energies of the 2^ states for Ca, 
Ti, Cr, Fe, Ni, Zn and Ge isotopes. The results except 
for Zn and Ge isotopes have already been discussed in 
our previous paper JTj. The lightest nucleus in each iso- 
tope chain is taken to be N—Z (cases with N < Z have 
mirror nuclei similar properties). The overall description 
of the 2+ energy levels is reasonable for all these isotope 
chains, although the calculated energies are systemati- 
cally higher than experimental ones by about 200 keV. 
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TABLE VI: Comparison of experimental electric quadrupole 
moments Qcxp (in unit efm^) with theoretical values Qth, 
which are calculated by using the effective charges = 1.5, 
= 0.5. Data are taken from Ref. |3^ . except for '' from 
Ref. 13] and which was evaluated based on * and a con- 
straint g(5*Fe; 10+)/O('^^Fe; 3/2-)=3.62±0.22 [H from the 
analysis of Mossbauer data. 



nuclei 


state 


(MeV) 


Qcxp 


Qth 


*''Ca 


7/2- 


O.QQQ 


+2.1(4) 


6.7 


*^Sc 


7/2- 


0.000 


-22(3) 


-20.6 




5/2" 


0.000 


+30.3(24) 


21.6 


«Ti 


2+ 


0.984 


-17.7(8) 


-12.6 


49^1 


7/2" 


0.000 


+24(1) 


22.0 


5«Ti 


2+ 


1.554 


+8(16) 


6.2 


50 y 


6+ 


0.000 


+20.9(40) 


19.0 


51 y 


7/2" 


0.000 


-4.3(5) 


-6.3 


50 Cr 


2+ 


0.783 


-36(7) 


-26.4 


52 Cr 


2+ 


1.434 


-8.2(16) 


-12.3 


53 Cr 


3/2" 


0.000 


-15(5) 


-15.3 


54 Cr 


2+ 


0.835 


-21(8) 


-24.4 


51 Mn 


5/2" 


0.000 


42(7) 


34.8 


52Mn 


6+ 


0.000 


+50(7) 


50.5 


54 Mn 


3+ 


0.000 


+33(3) 


33.2 


55Mn 


5/2" 


0.000 


+33(1) 


35.4 


5'»Fe 


2+ 


1.408 


-5(14) 


-22.6 




10+ 


6.527 


58(6)" 


53.5 


56Fe 


2+ 


0.847 


-19(8) 
16(1)'' 


-27.4 


5^Fe 


3/2" 


0.014 


15.8 


58Fe 


2+ 


0.811 


-27(5) 


-27.9 


5*^00 


4+ 


0.000 


+25(9) 


28.2 


5^Co 


7/2" 


0.000 


+52(9) 


36.0 




3/2" 


1.378 


+22(3) 


-22.5 




3/2" 


1.758 




13.8 


58 Co 


2+ 


0.000 


+22(3) 


22.3 


59 Co 


7/2" 


0.000 


+40(4) 


40.4 


«0Co 


5+ 


0.000 


+44(5) 


50.7 




2+ 


0.059 


+30(40) 


26.3 


58Ni 


2+ 


1.454 


-10(6) 


-2.4 




2+ 


1.333 


+3(5) 


3.9 


''iNi 


3/2" 


0.000 


+16.2(15) 


14.2 




5/2" 


0.067 


-20(3) 


-19.5 


62Ni 


2+ 


1.173 


+5(12) 


25.3 


64Ni 


2+ 


1.346 


+35(20) 


10.9 


•^^Cu 


3/2" 


0.000 


-21.1(4) 


-20.4 


''5CU 


3/2" 


0.000 


-19.5(4) 


-19.0 


«3Zn 


3/2" 


0.000 


+29(3) 


20.8 


i^^Zn 


2+ 


0.992 


-32(6) or -26(6) 


-7.4 


155 Zn 


5/2" 


0.000 


-2.3(2) 


-4.8 


6^Zn 


5/2" 


0.000 


15.0(15) 


16.1 


^"Zn 


2+ 


0.885 


-23.3(22) 


-2.5 


«^Ga 


3/2" 


0.000 


19.5 


19.6 


«»Ga 


1+ 


0.000 


2.77(14) 


-1.3 


e^Ga 


3/2" 


0.000 


+16.8 


17.5 


^iGa 


3/2" 


0.000 


+10.6 


12.8 


e^Ge 


5/2" 


0.000 


2.4(5) 


6.8 


^»Ge 


2+ 


1.039 


+3(6) 


19.8 


^^Ge 


2+ 


0.834 


-13(6) 


13.2 




Q(cal.) (efm") 

FIG. 12: Comparison of experimental electric quadrupole 
moments with the shell-model results. All data in Table IVII 
are included for which the sign is measured experimentally. 



In all cases, the energy jump corresponding to = 28 
shell closure is nicely reproduced for Ca to Ni isotopes. 
In the case of Zn isotopes, Ex{2f) is almost constant, 
which is also well reproduced. Recently measured value 
for ^^Ti (7V=32) comes precisely on the prediction of 
GXPFl. The new data of ^^Ct (iV=34) ^ also follows 
the predicted systematics. 

Very recent data for ^^Ti gives a 2+ energy of 1.13 MeV 
fisl l. This is significantly lower than the GPFXl pre- 
diction of 1.52 MeV. The KB3G interaction gives 0.89 
MeV for the '^^Ti 2+, and experiment lies in between 
GPFXl and KB3G. The dominant neutron component 
of this 2+ state has one neutron in the /5/2 orbit. In 
^^Ti there are some high spin states whose wave func- 
tions are dominated by the configuration with one neu- 
tron in the orbit, and their experimental energies 
are about 400 keV lower than GXPFl and in between 
GPFXl and KB3G The imphcation is that the ef- 
fective single-particle energy for the neutron /2 orbit in 
Z=22 is about 800 keV too high compared with GPFXl. 
The strong monopole interaction between the proton 
orbit and the neutron /5/2 orbit is responsible for 
lowering the energy of the /2 neutron energy (relative 
to P3/2) as protons are added to the /7/2 orbit from Z=20 
to 28 (see the right-hand side of Fig. 1 in 0). Thus to 
improve the agreement for ^^Ti, one would need to re- 
duce the strength of these two-body matrix elements in 
a way which is consistent with the entire fit. This will 
be one of the considerations for a next generation inter- 
action. The A'^=34 gap between pi^2 and f^/2 is about 
4 MeV. If it is estimated too large by 800 keV by the 
GXPFl interaction, the real gap turns out be greater 
than 3 MeV. Thus, this very new data seem to support 
the appearance of 7V=34 gap, while details of the GXPFl 
interaction may have to be improved. One must be also 
aware of the small separation energy of /5/2 in some of 
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the nuclei being discussed here. Such a small separa- 
tion energy induces additional relative lowering of /5/2, 
whereas this effect becomes weaker in nuclei with more 
protons. This effect is not included in the present fit, 
because it is not sizable in the nuclei used for the fit. 

The 2+ level for ^OCr (iV=36) ^ deviates toward 
lower energy as compared to the GXPFl prediction. This 
is correlated with the deviation in binding energy and 
both are signatures of the more collectivity from mixing 
with 59/2- In the Fe isotopes, a similar deviation can be 
seen at A^=38. These results suggest the limitation of the 
reliability of GXPFl interaction for neutron-rich nuclei. 
It is likely that one will need to explicitly introduce the 
gg/2 orbit into the model space in order to improve the 
calculations as one approaches N—AO. 

In the right panel of Fig. 1131 the E2 transition ma- 
trix elements B(E2;2^ 0^) are shown. Experimental 
B(E2) are significantly larger than theory for < 24, 
especially for the Ca isotopes, although the agreement in 
the excitation energy of 2^ state is reasonable. These de- 
viations show the large effects of the core excitation that 
has long been known for nuclei like ''^Ca 50]. The ■*°Ca 
core is significantly broken and we should take into ac- 
count the excitation from the sd-shell explicitly in order 
to reproduce these B(F,2) values. On the other hand, 
the agreement between theory and experiment is quite 
good in the middle of the shell, especially for Fe and 
Ni isotopes. The dependence on the neutron number is 
nicely reproduced including the A'^ = 28 magic number. 
For Zn and Ge isotopes with > 34, we again see the 
need of more collectivity in the model space, which also 
suggests the necessity of the 59/2 orbit, consistently with 
the results of electro-magnetic moments discussed in the 
previous subsections. 



F. Semi-magic nuclei 

In this subsection, spectroscopic properties are stud- 
ied in detail for several or Z = 28 nuclei around ^^Ni. 
If we assume an inert closed core for ^^Ni, these semi- 
magic nuclei are described by a few valence nucleons of 
one kind, i.e., only proton holes or only neutron parti- 
cles. Therefore, the low-lying level density is rather low 
and the effects of core-excitations are easy to interpret. 
In fact, it has been pointed out in Ref. that, al- 

though the shell-model calculation in a truncated t < 2 
subspace is quite successful for describing most of the 
low-lying stats in A^=28-30 nuclei, it is impossible to re- 
produce those states which contain sizable broken-core 
components, such as the lowest excited states in ^^Co, 
excited states above the lowest 3/2~, 5/2~, l/2~ triplets 
in ^''Ni, excited 0+ in ^**Fe and ^®Ni, etc. Similar difficul- 
ties can be seen also in Ref. [H^, where the shell- model 
calculations in a truncated space have been carried out 
with a different effective interaction. Note that the cause 
of this problem lies not only in the truncation of the 
model space but also in the effective interaction itself. 
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FIG. 13: (Left) First 2^ energy levels as a function of the 
neutron number A''. Experimental data are shown by filled 
circles square ?43|, open circle "4^, triangle (down) fisll . 
and triangle(up) [Slj. Solid lines show results of shell-model 



calculations. The truncation order t is 5 for 



56,58 



Fe, °"Ni, 



60,62,642jj^ 6 for 56'58.62Ni and 7 for ^'^•^^Fe. The other resuhs 
are exact. (Right) The values of _B(E2;_2j'" Of). Experi- 



mental data are shown by filled circles jSTjl . squares |39|], tri- 
angle (up) |53, triangles (down) [S^l, and open circles [i^. 



It is well known that the KB3 interaction (and also its 
descendants) gives an excellent description for low-lying 
states of A < 52 nuclei, but KB3 fails around ^^Ni, even 
when the model space is sufficiently large for convergence. 
Thus it is interesting to investigate whether GXPFl can 
properly describe these low-lying states which are sensi- 
tive to the core-excitation. We consider semi-magic nu- 
clei with Z=25, 26, 27 and A^=29, 30, 31. 
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FIG. 14: Experimental (right) and calculated (left) energy 
levels of ^''Mn. Experimental data are taken from [5fll |. The 
label "na" indicates that the spin-parity is not assigned exper- 
imentally. The calculated yrast states are connected with ex- 
perimental counterparts (or their candidates) by dotted lines. 



Figure El shows experimental and calculated energy 
levels of ^'^Mn. Most of the theoretical energy levels are 
obtained in the t = 7 subspace. The exact (t=13) results 
are also obtained in some cases. All experimental energy 
levels below 3 MeV excitation energy as well as the yrast 
states are shown with their theoretical counterparts and 
several additional states. 

For the yrast states, the agreement between the theory 
and the experiment is quite good. Assuming an inert '^^Ni 
core, this nucleus is described by three neutron holes in 
the f'j/2 orbit. Under such an assumption, possible states 
are J'' = 3/2", 5/2", 7/2", 9/2", 11/2" and 15/3". In 
the calculated wave functions, the lowest (/7/2)^^ con- 
figuration is in fact dominant (46^60%) for the yrast 
states of these spin-parities. Therefore, the truncation of 
the model space by a small t ~2 is reasonable for these 
states. 

Other yrast states, 1/2", 13/2", 17/2-, 19/2", 
21/2", 23/2" and 25/2" consist mainly of 
■^{hl2fv{hl2?{v?.l2t and T^{h I2fv{j-j I2)\!^I2Y 



configurations. These states are also well reproduced 
in the calculation, indicating that the effective single- 
particle gap above = 28 is reasonably reproduced in 
the present interaction. 

In order to investigate the properties of the core in de- 
tail, it is important to consider non-yrast states which 
contain a large amount of multi-particle excitations from 
the A^ = 28 core. Above E.,, ~ 2.2 MeV, the experi- 
mental level density increases significantly, which is well 
reproduced in the calculation. The largest difference be- 
tween the experiment and the calculation is found for 
the 3/22" state which lies at lower energy by 450 keV in 
the calculation, resulting in the inversion of 5/2^ and 
3/22". Nevertheless, one can find reasonable one-to-one 
correspondence between the experimental and theoretical 
energy levels. 

Focusing on the neutron configuration, the main con- 
figuration of these non-yrast states is v{h 12)^ (j'Y V 
stands for (j53/2,/5/2jPi/2)]- The typical probability for 
the neutron Ip-lh configuration is '^60%. These states 
contain also multi-neutron excitations as 2p-2h ('^25%) 
and 3p-3h (~7%). Therefore, taking into account the 
proton excitation, the i ~ 4 subspace is reasonable for 
describing these states. 

However, even in this energy region of Ex ^ 2.5 MeV, 
there are several states in which the A^ = 28 core is 
more severely broken. In 5/22^, Ijl^, 9/22^, the main 
neutron configuration is 2p-2h type (~60%), and also 
there are sizable broken-core components such as 3p-3h 
(~17%) and 4p-4h (~6%). Evidently, the < = 4 sub- 
space is no longer sufficient and at least t — % va needed 
to describe these states properly. In fact, for exam- 
ple, the electric quadrupole moment for 5/2^ varies as 
Q(5/2^)=-0.25, -1-0.32, -H0.36 and -h0.37eb for t=4, 5, 
6 and 13 (exact), respectively. The change of the sign 
in going from t = 4 to t = 5 is due to the crossing of 
5/2^ and 5/2^. A similar crossing is found in the case of 
7/2^. Although these neutron 2p-2h states are connected 
by relatively large (collective) Ml and E2 transition ma- 
trix elements as S(E2;7/2^ 5/2^) = 34 W.u. and 
B(Ml;7/22" 5/22") = 0.17 W.u., it may be difficult 
to observe such transitions because of negligibly small 
gamma-ray energies in comparison to the transition to 
the yrast states. 

In Table IVIII calculated electro-magnetic transition 
probabilities are compared with experimental data. Rea- 
sonable agreement can be seen in most cases, including 
transitions between non-yrast states. Note that the the- 
oretical 1/2^ is assigned to the experimental l/2~ at 
2671 keV so that the transition data can be consistently 
described. 
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In Fig. 1151 the energy levels of ^"^Fe are compared 
with experimental data. All experimental energy levels 
below Ex=4: MeV and all yrast states are shown with 
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TABLE VII: B(M1) and B(E2) for ^^'Mn. The excitation 
energyJEx) is shown in keV. Experimental data are taken 
from W\. 



initial 


final 


multi- 


exp. 


th. 




J"(£.) 


pole 


(W.u) 


(W.u) 


5/2" (378) 


7/2- (0) 


Ml 


0.00254(23) 


0.0060 






E2 


14(3) 


21.7 


3/2^(1290) 


5/2- (378) 


Ml 


0.023(3) 


0.0194 






E2 


1.9(5) 


0.6 




7/2- (0) 


E2 


13.4(11) 


8.1 


11/2" (1441) 


7/2- (0) 


E2 


12.8(18) 


9.4 


9/2" (1620) 


5/2- (378) 


E2 


3.7(6) 


3.6 




7/2- (0) 


Ml 


0.0012(3) 


0.0002 






E2 


7.0(9) 


5.4 


5/2" (2274) 


3/2- (1290) 


Ml 


0.0021(8) 


0.0000 






E2 


3.0(15) 


0.0 




5/2- (378) 


Ml 


0.0019(4) 


0.0010 






E2 


0.61(14) 


1.6 




7/2- (0) 


Ml 


0.0054(11) 


0.0023 






E2 


0.07(4) 


0.02 


3/2" (2407) 


3/2- (1290) 


Ml 


0.064(24) 


0.0482 






E2 


1.9(16) 


0.3 




5/2- (378) 


Ml 


> 0.0038lio 


0.0058 






E2 


< 0.43 


1.0 




7/2- (0) 


E2 


2.0(9) 


4.4 


13/2- (2563) 


11/2- (1441) 


Ml 


0.00146(18) 


0.0022 






E2 




0.0 


1/2" (2671)" 


5/2- (378) 


E2 


15(12) 


2.7 


7/2^(2686) 


5/2- (378) 


Ml 


0.004(4) 


0.0014 






E2 




0.1 




7/2- (0) 


Ml 


0.010(9) 


0.0100 






E2 


2.4(22) 


0.9 


15/2" (2693) 


13/2- (2563) 


Ml 


0.21(4) 


0.0183 






E2 




4.7 




11/2- (1441) 


E2 


5.4(8) 


3.3 


3/2" (2876) 


5/2- (378) 


Ml 


2.5E-5(10) 


0.0040 






E2 




0.9 




7/2- (0) 


E2 


0.0008(3) 


0.2 


9/2^(2947) 


7/2- (0) 


Ml 


0.007(4) 


0.0315 






E2 


1.8(8) 


3.9 




1 1 /9'" ^'^«Q9i^ 

ll/Z yZOiJoj 


Ml 


n 0/1 

xJ.yJ^yo ) 








E2 


1.6E+2(12) 


18.1 


15/2- (3439) 


15/2- (2693) 


Ml 


0.23(4) 


1.0585 






E2 


2.E+1(5) 


1.3 


19/2- (5614) 


17/2- (4384) 


Ml 


> 0.22 


0.1183 






E2 




16.1 


21/2- (6533) 


17/2- (4384) 


E2 


2.8(16) 


7.0 


23/2- (7004) 


21/2- (6533) 


Ml 


0.153(24) 


0.1099 






E2 




0.0 
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their theoretical counterparts. Shell-model calculations 
have been carried out in the t — 7 or larger subspace 
including the exact (t=14). One can see remarkable one- 
to-one correspondence between theory and experiment 
for most of these levels. 

The yrast band 0+-2+-4+-6+ shows a typical proton 
two hole spectrum in the /7/2 orbit under the influence 
of a short range interaction. In fact the wave functions of 
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FIG. 15: Energ 
from 1^ and y 



levels of ^''Fe. Experimental data are taken 
|. Conventions are the same as in Fig. 1141 



these states are dominated by the TT{f7/2)^'^{h/2)^ con- 
figuration (50-60%). On the other hand, the yrast 8+ and 
10+ states consist of core-excited configurations. While 
the dominating configurations are still of neutron Ip-lh 
excitation type from the N = 28 core (36% and 57% for 
8+ and 10+ states, respectively), the wave functions can 
no longer be approximated by only one configuration. 
The energy gap between 6+ and 8+ reflects the softness 
of the ^^Ni core, which is determined by the effective in- 
teraction. It is shown in Ref. "56] that the energy gap 
is too large with the KB3 interaction. This is consistent 
with the fact that it also gives an energy for the yrast 2+ 
energy in ^^Ni which is too high compared with experi- 
ment. 

It was pointed out in Ref. that the FPD6 

interaction predicts the existence of deformed states at 
relatively low excitation energy (^3 MeV), which consist 
of neutron 2p-2h configurations. In the present calcula- 
tion, the candidates of such states appear as O2 and 23". 
The neutron 2p-2h configuration is 35% and 26% in these 
states, respectively, which are the most dominant compo- 
nents in their wave functions. However, the deformation 
of these states is much smaller than the FPD6 prediction. 
The electric quadrupole moment for the "2p-2h" 2+ state 
is predicted to be — 0.30eb by the variation after the 
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TABLE VIII: B{M1) and B(E2) for ^''Fe. Experimental data 
are taken from j57|. 



initial 


final 


multi- 


exp. 


th. 
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sition matrix elements are compared with experimental 
data. The E2 transition matrix elements in the ground- 
state band are reproduced fairly well in the present cal- 
culations, although the 8+ 6+ transition is overesti- 
mated. Such a deviation is commonly seen also in the 
results of KB3 and FPD6 interactions 1>%. Other B(E2) 
values are reasonably reproduced except for several of 
the small ones. As for the i3(Ml) values, the agreement 
between theory and experiment is basically good. The 
exception is 4+ (3294) 4+ (2538) transition, which is 
significantly overestimated by more than one order of 
magnitude in the present calculation. 



angular momentum projection (VAP) method with the 
FPD6 interaction. On the other hand, in the present 
results, it is (5(23") — —0.09 eb, and the E2 transition 
matrix element is B{E2; 2^ Oj) = 11.3 W.u., which is 
slightly larger than that of the yrast in-band transition. 
These values are inconsistent with the axially symmet- 
ric large prolate deformation. In the present results, two 
neutrons are excited mainly to the orbit, while it is 
/5/2 in the results of the FPD6 interaction, reflecting the 
fact that the effective single particle energy for the 
orbit is too low for FPD6. 

In Table IVIIII the calculated electro-magnetic tran- 



3. Co 

In Fig. ^1 energy levels of ^^Co are shown. All ex- 
perimental energy levels below Ex=^ MeV and all yrast 
states up to J'^ = 23/2 ~ are compared with theoretical 
results. There is a state with unknown spin-parity at 
i?2;=2.960 MeV, which can be interpreted tentatively as 
a 5/2^ state in the correspondence to the present calcu- 
lation. The overall agreement between theory and exper- 
iment is quite good. 

If we assume an inert ^^Ni core, this nucleus is de- 
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scribed as one neutron hole in the /7/2 orbit, correspond- 
ing to only one 7/2~ state. In the calculated ground-state 
wave function, the probability of this Op-lh configuration 
relative to the ^^Ni core is only 65%, and there are sizable 
(^20%) 2p-3h configurations. All excited states consist 
of core-broken configurations even in the lowest order of 
the approximation. Therefore the property of the core is 
expected to manifest clearly in the yrast spectrum. 

The yrast states up to J'^ = 17/2^ except for the 
ground state are described mainly by lp-2h configura- 
tions (30~56%). Thus the truncated subspace with small 
t ^ 2 is expected to be reasonable. In fact the truncated 
t = 2 subspace was used in the shell-model calculations 
of Ref. [SJI, where 13/2~ and higher spin states are de- 
scribed reasonably well. However, the excitation energies 
of lower spin states 9/2~ and ll/2~ are predicted to be 
too high by ~1 MeV, and the results of 5/2~, 3/2~ and 
l/2~ are not shown in Ref. 

A part of such a difficulty may be attributed to the 
severe truncation. In the present t — 7 results, the po- 
sition of both these higher (J'^ > 13/2~) and lower spin 
states are successfully described. In the lower spin states, 
there are sizable 3p-4h (^^21%) and 4p-5h (~9%) com- 
ponents, including both proton and neutron excitations 
almost equally. On the other hand, the higher spin states 
13/2~, 15/2^ and 17/2~ are described mainly by neutron 
excitations only, and the probabilities of 3p-4h and 4p- 
5h components are much smaller (~16% and 6%, respec- 
tively). When the GXPFl interaction is used in the t = 2 
subspace, the excitation energies of lower spin states be- 
comes too large by about 0.4 MeV, while those of higher 
spin states are almost unchanged. 

For constructing higher spin states, it is necessary to 
excite one more nucleon from the /7/2 orbit. In fact the 
most dominant component in the wave function of these 
states is the 2p-3h type ('-45%). A large gap (~2 MeV) 
between 17/2" and 19/2" reflects this core-excitation. 
Note that, although it is possible to construct a 19/2~ 
state within the lp-2h configuration space, this compo- 
nent is only 9% in the 19/2 state, similarly to the result 
in Ref. (5^]. The present calculations successfully repro- 
duce both lp-2h and 2p-3h excitations. 

The structure of non-yrast states is more complicated. 
Since the present interest is the low-lying core-excited 
states, we can focus on 1/22^, 3/23", 7/23 and 5/2^. 
These states are all below -Ea;=4 MeV and contain more 
than 45% neutron 2p-2h components. A similar charac- 
ter can also be seen in 3/22^, although the probability of 
the neutron 2p-2h component is slightly smaller (37%). 
These states are connected by relatively large E2 transi- 
tion matrix elements such as i3(E2; 7/2;^ — > 3/23" )=11 
W.u. and B{E2; 5/2^ ^ 1/2^)=18 W.u. It wiU be a 
good test if it is possible to observe such transitions. 

The calculated 3/24 has isospin T = 3/2 and is the iso- 
baric analog state of the ^^Fe ground state. The correct 
position of this state, as seen in comparison to experi- 
ment in Fig. 1161 implies the proper isospin structure of 
the present interaction. 



TABLE IX: B{M1) and B(E2) for ^^Co. Experimental data 
are taken from 58]. 
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"5/23 in the calculation. 



In Table IIXI electro-magnetic transition matrix ele- 
ments are shown for ^^Co. We can find reasonable agree- 
ment between theory and experiment. In general, the 
deviations are large where the experimental errors are 
also large. 

4. ^^Ni 

In Fig. El calculated energy levels of ^^Ni are com- 
pared with experimental data. All calculated states up 
to Ex=& MeV as well as the yrast states are shown. The 
shell-model results were obtained in the t — 7 subspace. 
The agreement between the experiment and the calcula- 
tion is basically good especially for even spin yrast states. 
However, the yrast 3+ and 5"'" states are not found in the 
experimental data, while no candidate for the experimen- 
tal (2+), (4+) and (6"'") can be seen in the calculation at 
least around ~6 MeV. More precise experimental in- 
formation is needed to discuss such discrepancies. Similar 
results are shown in Fig. 2 of Ref. for even spin states, 
which were obtained by the MCSM calculations, taking 
typically 13 bases per one eigenstate. We thus confirm 
the reliability of the MCSM results for non-yrast states. 

For the yrast states, the ground state consists mainly of 
the (77/2)^® closed-shell configuration (68%). The excited 
states 2+, 3"*", 4+, S"*" and 6+ are of the Ip-lh character 
(42^50%) relative to the closed-shell configuration, and 
the excitation energy of 2^^ reflects the shell gap at N or 
Z — 28. The 2p-2h conflgurations are dominant in 1+, 
7+, 8+, 9+ and 10"'" states (35~51%). One can flnd again 
an energy gap between these Ip-lh and 2p-2h groups. 
The 3p-3h dominant states appear above ~10 MeV. 
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FIG. 18: Convergence of the calculated energies (lower panel) 
and quadrupole moments (upper panel) as functions of the 
truncation order t for the lowest two 2"^ states of ^''Ni. Hori- 
zontal lines show the corresponding MCSM results. 



In the same way, most of the non-yrast excited states 
shown in the figure can be classified as either Ip-lh (3^, 
5^, • • • ) or 2p-2h (2^, 6^, • • • ) dominant states. These 
np-nh states contain typically 20% of (n + 2)p-(n -I- 2)h 
configurations. Therefore, the 4p-4h configuration space 
is essential (but not enough) to describe these states. 

As has already been discussed in Ref. 0, GXPFl pre- 
dicts the deformed 4p-4h band as O;^, 2^, • ■ • . In the 
present results, 0;^ consists mainly of 4p-4h (48%), 5p-5h 
(23%) and 6p-6h (12%) configurations. The structure of 
2j is quite similar. Special care should be taken con- 
cerning the convergence of the calculation for such a de- 
formed 4p-4h band, since the present truncation scheme 
is obviously not suitable for describing strongly deformed 
states, as discussed in Ref. IlSj for the case of the FPD6 
interaction. Figure El shows the calculated energies and 
quadrupole moments of the lowest two 2+ states as a 
function of the truncation order t. It can be seen that 
the 4p-4h state 2^ is not completely converged even in 
the t = 7 subspace, which corresponds to the M-scheme 
dimension of 89285264. In the same figure, the MCSM 
results are also shown, where 25 bases were searched for 
each eigenstate. The MCSM energy of the yrast 2+ state 
is almost the same as that of the t = 7 calculation, which 
shows reasonable convergence. On the other hand, for 2j 



states, the MCSM energy is lower than the t = 7 result 
by about 0.2 MeV, and should be more accurate. The 
electric quadrupole moment is (3(2^) = — 41.2efm^ in 
the MCSM calculation, while it is —33.6 in the present 
t — 7 result. Similar results are obtained for the 4p- 
4h 0;^ state. The Af-scheme dimension of t = 8 subspace 
reaches to 255478309, which is beyond our computational 
limitations. For the description of such deformed states 
where a lot of mixing among various spherical config- 
urations occurs, the MCSM calculation is more efficient 
and suitable, because certain basis states can be sampled 
from and around local minima. 

In the present results, we can find several T = 1 states 
such as Oj, 24" and 4^. These states are interpreted as 
the isobaric analog states of ^^Co. The excited O"*" with 
T = 1 is in fact found experimentally at excitation energy 
of 7.904 MeV, which is consistent with the present result. 

Experimental data for the electro-magnetic transitions 
are rather limited. We can only compare with the few 
transition strengths shown in Table Ixl 
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TABLE X: B(M1) and B(E2) for ^''Ni. 
are taken from and j52l] . 
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FIG. 19: Energy levels of ^^Ni. Experimental data are taken 
from |6l| and |54|. Conventions are the same as in Fig. 1141 



Energy levels of ^^Ni are shown in Fig. ^| All exper- 
imental yrast levels as well as non-yrast states up to 4 
MeV excitation energy are compared with the theoreti- 
cal predictions. The shell-model calculations were carried 
out in the t = 6 subspace. The agreement between ex- 
periment and theory is basically good up to high spin 
states 23/2~, although there are several inversion of the 
order of two closely lying states, such as 7/2 j" and 5/22^. 

Comparing the present results with the previous 
MCSM ones in Ref. 0, we can again confirm the re- 



liability of the MCSM calculations for non-yrast states 
as well as the yrast states of odd-mass nuclei. The root 
mean square difference of the excitation energies between 
the present and the previous MCSM results is 154 keV 
for 19 states shown in Fig. 2 of Ref. 0. 

If we assume an inert ^^Ni core, ^^Ni is described by 
only one neutron. Then, three states 3/2~, 5/2" and 
l/2~ are possible, corresponding to the occupation of 
P3/2, /5/21 a-iid Pi/2 orbits by one neutron. According 
to the present calculations, the lowest three states show 
such a "single particle" character, although there is a 
sizable mixing of the broken-core configurations. In fact, 
the pure "single particle" configurations are 62, 58 and 
46%, respectively, in the calculated wave functions. Since 
the closed core configuration in the ground state of ^^Ni 
is 69% 7] , the core is further broken in the ground state 
of ^^Ni, as can also be seen in Fig. |51 

The bare single particle energy of the /5/2 orbit rela- 
tive to the P3/2 orbit is found to be 4.3 MeV in GXPFl. 
(see TablcQJ. The effective single particle energy (ESPE) 
[5IITH of the /5/2 orbit relative to the P3/2 orbit decreases 
rapidly as the proton 77/2 orbit is occupied, and reaches 
1.1 MeV at ^^Ni. The mechanism for this behavior of 
the orbit is explained by the strong attraction be- 
tween the protons in the fj/2 orbit and neutrons in the 
/5/2 orbit due to a strong j> - j< coupling term in the 
nuclear force, as discussed in Ref. 0. Thus the proper 
contribution of the two-body interaction, especially the 
monopole part, is needed to reproduce the experimental 
order of the lowest three states in ^^Ni. The FPD6 in- 
teraction is too strong in this respect, giving an ESPE 
for the /5/2 orbit which is too low for the Ni isotopes. 
Note that, the ESPE's of the /5/2 and the pi/2 orbit are 
still too large by 0.3 and 0.9 MeV, respectively, in com- 
parison to the experimental excitation energies of 5/2~ 
and 1/2" states, since they are defined in terms of a 
closed-shell structure for ^^Ni. Further energy gain by 
the mixing of the broken-core components is needed for 
the better reproduction of the experimental data. 

There is an energy gap (~ 1.3 MeV) above the low- 
est three states, which corresponds to the 2p-lh core- 
excitation. The doublets 5/22^ and 7/2j~ consist of two 
dominant configurations, i^ifr /lY {Pz/2)^v{f f7 /2f{P3/2)'' 
and 7r(/7/2)*J^(//7/2)''(p3/2)^, sharing almost equal prob- 
abilities (^20%). Above these doublets, most of the 
states up to ~5 MeV are also of the 2p-lh nature, 
although one can no longer find single dominant config- 
uration which exhausts more than 20% probability. In 
this 2p-lh regime, as exceptional cases, the 3/23" and 
1/22^ states are comprised mainly of higher configura- 
tions: the 3p-2h configuration takes the largest weight 
(25 and 29%, respectively) and even the 5p-4h configura- 
tion can be seen with a non- negligible probability ('^5%). 
In the yrast 15/2" and most of the higher spin states, the 
3p-2h configuration becomes dominant, carrying more 
than 40% probability. We have not identified 4p-3h or 
more significantly core-broken states such as those of the 
deformed 4p-4h band in the ^^Ni in the present results 
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TABLE XL B(M1) and B(E2) for "Ni. 
are taken from Ifill. 
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for 57Ni. 

Electro-magnetic transition strengths are shown in Ta- 
ble IXII Because of the large ambiguities in the experi- 
mental data, it is difficult to draw a definite conclusion. 
Most of the calculated Ml transitions are consistent with 
the experimental data. A notable deviation is found in 
the B{M1; 5/2j~ — > 3/2j~), where the calculated value is 
too small by a factor of 8. As for the E2 transitions, it can 
be seen that the calculated i3(E2) values are in general 
smaller than the experimental data. The exception is the 
transition from 5/2- (2443). However, no experimental 
error is indicated in the data and only the upper bound 
is given for the associated Ml transition. The calculated 
i3(E2) value does not contradict with the experimental 
lifetime data. 



Figure 1201 shows calculated and experimental energy 
levels of ^^Ni. The yrast states up to 10 MeV excitation 
energy, non-yrast states below Ex=^ MeV, and several 
additional states are shown with their theoretical coun- 
terparts. The agreement between the experiment and the 
calculation is satisfactory. The calculations were carried 
out in the t=6 truncated subspace. The results are ba- 
sically consistent with our previous MCSM calculations 
in Ref . jTj , although the latter spectrum was slightly ex- 
panded for higher spin states due to the small number of 
basis states (~13 per one eigenstate). 

In the calculated ground-state wave function, the 
total probability of the closed core configurations 
''■(/7/2)^f^(/7/2)^(^)^ is only 57%, which is smaller than 
those in ^®'^^Ni. As shown in Fig. O this quantity de- 
creases further for larger N and takes a minimum value at 
= 34 (52%). Thus the ^^Ni core is soft in Ni isotopes. 
Nevertheless, for understanding the basic structure, it is 
still useful to assume an inert ^^Ni core and to consider 
two neutrons in the P3/2, /5/2 f^nd pi/2 orbits. Then, the 
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FIG. 20: Energy levels of ^®Ni. Experimental data are taken 
from Conventions are the same as in Fig. 1141 



allowed maximum spin is J = 4, and core-excitation is 
needed to generate higher spin states. In the calculated 
wave function, the leading configuration is v{p'i/2Y for Qi 
and 2f (37% and 30%, respectively), and v{p^/2Y{fb/2Y 
for 1^ and 3^ (42% and 52%, respectively). There is 
a large energy gap (~ 2 MeV) in the yrast spectrum 
between 4^ and 5^, reflecting the core-excitation. One- 
proton core-excitation dominates the wave functions up 
to J = 9, although there is no single configuration which 
exhausts more than 20% weight in the calculated wave 
functions. For J = 10 and higher spin yrast states one- 
proton one-neutron core-excitation becomes the most im- 
portant mode. Consistently, there is another energy gap 
(~ 1.6 MeV) above %i in both calculated and experimen- 
tal energy spectra. 

The closed-core configurations also play a major role 
in many non-yrast states. Among the non-yrast states 
below Ex ^ ^ MeV, the 3p-lh configurations associated 
with the one-proton core-excitation become dominant in 
2;^, •^2- Since the level density increases signifi- 

cantly above ^ 4 MeV, it is difficult to make a simple 
interpretation for the properties of the core-excitations 
in higher states. 

As discussed in Ref. 0, the core-excitation of two- 
protons is dominant in the wave function of Oj|". In the 
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TABLE XII: B(M1) and B(E2) for ^**Ni. Experimental data 
are taken from j62l]. 
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present calculation, we found that 2g shows a similar 4p- 
2h character. The probability of such a 4p-2h configura- 
tion Tr{h/2f{rfj^{h/2f{rf is 46 and 41% for 0;^ and 
2g", respectively. The calculated quadrupole moment is 
Q(2+) = -0.30 eb, and B(E2; 2+ ^ 0;[)=18 W.u. These 
values are consistent with those of the rotational model, 
corresponding to the K ~ Q band on top of the axially 
symmetric prolate rotor with an intrinsic quadrupole mo- 
ment Qo=\\{) cb. The deformation of the present 4p-2h 
band appears to be slightly smaller than the 4p-4h band 
in ^^Ni 7], which is consistent with Qo=140 eb obtained 
by using the same effective charges. 

In Table IXIll calculated transition strengths are com- 
pared with experimental data. In general, the calculated 
Ml matrix elements are too small. The exception is 
_B(M1; 2^ — > 2^), which is too large by one order of 
magnitude. As for the E2 matrix elements, the calcu- 
lated values are almost consistent with experimental data 
in the case of AJ = 2 transitions, i.e., no mixing of Ml 
transitions. One notable exception is i3(E2; Oj" — > 2]^), 
which is too large by more than 3 orders of magnitude, 
but this is an extremely weak transition. 



7. ^^Ni 

Energy levels of ^^Ni are shown in Fig. |^ All exper- 
imental data for yrast states as well as non-yrast states 
below 3 MeV excitation energy are compared with the 
shell-model results obtained in the t=Q subspace. It is 
found that the theoretical results agree quite well with 
the experimental data. The results of the lowest triplet 
3/2^, 5/2" and l/2~ have already been reported in Ref. 
7] . The present calculation has confirmed the validity of 
GXPFl interaction also for higher excited states. 

Assuming an inert ^^Ni core, this nucleus is described 
by three valence neutrons in the upper three orbits P3/2, 
/5/2 and Pi/2- Therefore various states are possible with- 
out any core-excitation. In fact such closed-shell (3p-0h) 
configurations carry the maximum weight (31%~56%) 
in most of the calculated low- lying states below E^— i 
MeV. There are three exceptions, 7/2^, 7/2;^ and ll/2j^, 
where 4p-lh configurations take the largest weight. Such 
4p-lh configurations also dominate the wave functions of 
higher-lying excited states shown in the same figure. The 
5p-2h configurations appears to be the leading configu- 
rations in the yrast 19/2" and 21/2" states. 

Like other nuclei, we would like to identify states with 
large components of 5p-2h or higher configurations at 
Ex <^ i MeV. This energy boundary is set because the 
level density increases very fast above it in ^^Ni. How- 
ever, as seen from the above discussions, there are no such 
candidate states theoretically or experimentally. There- 
fore, the study of 5p-2h and higher states with lower an- 
gular momenta is not possible in ^^Ni. Thus, the dis- 
cussion of the multi-particle multi-hole states which is 
one of the main elements in our assessment of the effec- 
tive interaction is confined mainly to a few nuclei around 
^^Ni. 

Table IXTTTl shows the measured and calculated B(E2) 
and B{M1). Even below Ex=2 MeV, there are several 
uncertain levels in the experimental data at 1695 keV 
(no spin assignment), 1739 keV (assigned to (9/2~)), and 
1746 keV (assigned to 5/2", 7/2"). In the present calcu- 
lations, we have tentatively assigned no theoretical states 
to these three uncertain states, since no probable candi- 
dates appear around E^ ^ 2 MeV. More experimental 
information is needed for a detailed comparison between 
theory and experiment. 

In general, the order of magnitude is well reproduced 
for most transitions. More specifically, the calculated 
Ml transition matrix elements are small in compari- 
son to the experimental data in many cases, except for 
1/2]" 3/2]" transition. The calculated E2 matrix ele- 
ments are also smaller than the experimental data. How- 
ever, one can clearly see the correlation between the cal- 
culated values and the corresponding experimental data. 
The agreement could be improved if we introduced larger 
effective charges. 
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FIG. 21: Energy levels of ^^Ni. Experimental data are taken 
from |6^ . Conventions are the same as in Fig. 1141 



TABLE XIII: B(M1) and B(E2) for ^''Ni. Experimental data 
are taken from 
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G. N=Z odd-odd nuclei 

In the p/ shell, N = Z odd-odd nuclei have been of 
special interest, where the lowest isospin T=0 and T = 1 
states are almost degenerate near the ground state. In 
most cases T 1, J = is the ground state. The only 
known exception is ^^Cu, in which the ground state is 
T = 0, J = 1. This reflects a detailed interplay of T = 
and T = 1 interactions, and it is important to evaluate 
the effective interaction from this viewpoint. 



1. '^^Co 

In Fig. 123 calculated energy levels of ^"^Co are com- 
pared with experimental data. The spin-parity has not 
been assigned to most of the states above Ex — 3 MeV, 
and there are several uncertain levels among low-lying 
states. It can be seen that the present shell-model calcu- 
lations in the t = 6 subspace give a reasonable descrip- 
tion, although there are several small deviations such as 
the interchange of the order of 2^ and 3^. The calcu- 
lated 0^2 J '^t^ "^3.4 ^'^^ ^2" T — 1, corresponding to 
the isobaric analog states of ^^Fe (see Fig. I15|l . while 
other calculated states shown in the figure are T = 0. 



In the calculated wave functions of the yrast odd-spin 
states 1^, 3^, 5^ and 7^, the lowest t^{!-j j-if v^j'j 
configuration is dominant (59, 58, 39 and 62%, respec- 
tively) , but core-excited components are not small. Such 
lowest 0p-2h configurations relative to ^^Ni core are also 
dominant in even-spin 0]^, 2+, 4+, and 6^, which are all 
T = 1. 

In the shell-model calculations using the Surface 
Delta interaction in the (/7/2, P3/2) space with a restric- 
tion that the P3/2 orbit can be occupied by one proton 
and one neutron, the 4+ and 6+ states of T = 1 appear as 
yrast states. On the other hand, 4^2 ^^re isospin T = in 
the present calculation. If the calculated 4^ is assigned 
to the experimental state (3"^, 4+) at 2083 keV, the agree- 
ment between theory and experiment is good. This state 
consists mainly of the lp-3h configurations (50%). A 
similar structure can be seen for 6^, 5^, 7^, and ^X^i 
where one nucleon is excited across the shell gap mainly 
to the P3/2 orbit. The structure of 4^ is similarly of the 
lp-3h type, although the excitation to the /5/2 orbit is 
dominant. Such an excitation is also large in 8^. 

One can find a state in the calculated spec- 
tra at a reasonable excitation energy, which does 
not appear in the shell- model results of Rcf. [6^ . 
The most important configuration in this state is 
l2f{V3l2fv{j'j i2f{j>zliY , but its probability is only 
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TABLE XIV: B(M1) and B(E2) for ^^Co. 
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15%. Such 2p-4h configurations are 44% in total, and 
there are other sizable core-excited components such as 
the 3p-5h (32%) and 4p-6h (17%) types. Similar core- 
excitations also appear in 2^, 3 J and 0^ (T = 1). 

The calculated Ml and E2 transition strengths are 
listed in Table IXTVl Since no experimental value is avail- 
able, the present results are compared with those of the 
shell- model calculation in Ref. 64]. Both results agree 
within a factor of about 2 for transitions between the 
states which consist mainly of the 0p-2h configurations. 
The strong isovector Ml transitions predicted based on 
the "quasideuteron" picture 0| are found also in the 
present large scale calculations. Note that the effective 
spin g-factors g'f^ — O.lg^g'^'^ were used in Ref. If 
we adopt the same g-factors in the present calculations, 
the i3(Ml) values are significantly reduced, for example, 
B(M1; 1+ 0+) becomes 1.25 W.u. Thus the configu- 
ration mixing strongly affects the _B(M1) value. 

Several resuhs related to the 1+ (1614), 3+ (2174) and 
4+ (2852) states are very different between two calcula- 
tions, because the structure of the corresponding states 
are different as mentioned above. 



2. 



'•Cu 



In Fig. 1221 the calculated energy levels of ^^Cu are 
compared with experimental data. All experimental 
states up to 2.2 MeV as well as the yrast states are 
shown. The shell-model calculations have been carried 
out in the t — & truncated subspace. The agreement be- 
tween theory and experiment is satisfactory. In the cal- 
culated spectrum, 0^, 2 J, 4^ and 0^ are isospin T — \. 
The Tz = 1 members of the multiplets can in fact be 
seen in the calculated spectrum of ^^Ni (see Fig. I2(J|I . 



4o in the calculation. 



The T — ground state 1^ is successfully reproduced. 

Assuming an inert ^^Ni core, one proton and one neu- 
tron are valence particles in the P3/2, /5/2 and Pi/2 or- 
bits. Therefore J < 5 states are possible. In the cal- 
culated wave functions, the yrast states up to 5"*" are in 
fact dominated by such closed core configurations with 
probabilities of 44 63%. The 6^ state consists mainly 
of 3p-lh configurations (45%). The 4p-2h configurations 
are dominant in 7^", 8^, 9f and IQ^ (41^59%), and the 
5p-3h configurations become the major components in 
llj*" and 12^ ('-40%). Thus it is natural that the t = 2 
shell-model results in Ref. |3 show apparent discrepan- 
cies from experiment for 11+ and 12+. 

In the present calculation, all states with J < 10 
are dominated by either 3p-lh or 4p-2h configurations. 
Therefore, shell-model calculations even in the severely 
truncated {t ^ 2) subspace give reasonable results at 
least near the yrast line. In order to obtain further in- 
formation on the core-excitation, it is necessary to study 
higher spin states or highly excited states far from the 
yrast line. From this viewpoint, recent precise measure- 
ments of the Gamow- Teller strength distribution in 
this nucleus up to high excitation energy provide quite 
important information. The GXPFl reproduces the mea- 
sured strengths around ~4 MeV reasonably well, 
which are missing in the prediction of the KB3G inter- 
action. 

In the calculated wave function of the ground state 
1]^, the configuration 7r(p3/2)^ 1^(^3/2)^ takes the largest 
weight (18%) and the second largest components are 
7r(P3/2)^J^(Pi/2)^ and 7r(pi/2) ^(^3/2)^ (13% for each). 
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TABLE XV: B(M1) and B(E2) for ^*Cu. Experimental data 
are taken from ^ a.nd^. 
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This is not very far from the "quasideuteron" picture 
although the configuration 7r(p3/2)^i^(P3/2)^ is dis- 
tributed among several 1+ states such as (15%) 
and l'^ (20%). However, it contradicts with the shell- 
model results in Ref. !54|, where 7r(p3/2)^i^(/5/2)^ and 
^(/5/2)^'^(P3/2)^ are dominant (20% for each). 

The calculated B{M1) and B(E2) values are compared 
with experimental data in Table IXVl Although the ex- 
perimental data are limited, the agreement between the- 
ory and experiment is reasonable. More detailed discus- 
sions are found in Ref. 1691. 



H. GXPFl vs. GXPFM 

In the previous subsections, it has been shown that 
GXPFl properly describes the low-lying core-excitations. 
As discussed above, the effects of the core-excitations can 
be observed directly in the low-lying energy spectra of 
semi-magic nuclei, and therefore the description of such 
core-excited states plays a crucial role for the evalua- 
tion of the effective interaction. From this viewpoint, in 
this subsection, we discuss a specific problem with the 
GXPFM interaction: an interaction derived by adding 
empirical corrections to the microscopic G interaction 
only in the monopole part and pairing matrix elements 
(see subsection III F|) . With the FDA* estimate, the de- 
scription of energy data by GXPFM appears to be poor 
especially for the yrare states of semi-magic nuclei. 

As examples, we consider the low-lying energy levels 
of s^Mn and ^"^Fe, which are shown in Fig. |^for yrast 
and yrare states. Details of the shell-model calculations 
are the same as those described in subsection IIII Fl For 
a comparison, the results of GXPFl are also shown. It 
can be clearly seen that, in both nuclei, the yrast states 
are described nearly equally well by both GXPFl and 
GXPFM at least in the region of < 3.5 MeV. On the 
other hand, several yrare states such as 5/2^ and 7/2^ 
in '"^^Mn and Oj in ^'^Fe are predicted to be too low by 
GXPFM. According to the analysis of the shell-model 
wave functions, these states are dominated by neutron 
core-excited configurations. Typically 40% of these wave 
functions consist of neutron 2p-2h configurations relative 
to the = 28 core. In the case of GXPFl, such states are 
predicted with a similar structure, while their excitation 
energies are closer to the experiment. 

This fact suggests that the "monopole + pairing" cor- 
rection works well for the description of the states in 
the lowest configuration and those with Ip-lh excita- 
tions relative to it, but it is not sufficient for treating 
the states which are dominated by the 2p-2h or more 
excited configurations. This is natural because, in the 
former states, the structure is expected to be "single- 
particle" like, especially for the semi-magic nuclei, and 
therefore such states can be described well by adjusting 
the monopole part only. On the other hand, in the latter 
states, the contribution of the multipole part should be 
more important, and multipole parts must be modified 
for their proper descriptions. Thus, the core-excitations 
provides us with a clue to investigate corrections to the 
interaction that go beyond monopole plus pairing. 



I. Quadrupole corrections 



In subsection III Bl it has been shown that major mod- 
ifications of the microscopic interaction can be found in 
the diagonal parts which are not necessarily of monopole 
character. We have also shown in the previous subsection 
that the "monopole + pairing" correction is insufficient 
for describing the property of core-excitations. Such in- 
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FIG. 24: Comparison of low-lying energy levels of (upper) 
^^Mn and (lower) ^'*Fe between the experimental data and 
the shell-model predictions by GXPFM and GXPFl. 



sufficiency can be seen explicitly in the energy spectra of 
semi-magic nuclei. 

Figure shows a part of the low-lying energy levels of 
56Ni and "Ni, where the results of GXPFl, GXPFM and 
KB3G are compared. The latter two interactions predict 
too high excitation energies for the 2f state in ^^Ni and 
the 5/2^ state in ^^Ni. Both of these excited states are 
dominated by the configuration with one nucleon excited 
from the /7/2 orbit to the P3/2 orbit. 

We have found that in this case the most relevant ma- 
trix elements are 1^(7373; Jl), which are compared in 
Fig. for G, GXPFl, GXPFM and KB3G. It can 
be seen that the difference between G and GXPFM is 
primarily of monopole character (constant shift). Simi- 
lar monopole corrections are present in KB3G, although 
its original microscopic interaction is not G but KB, re- 
flecting that G and KB are very close for these matrix 
elements. On the other hand, GXPFl shows a rather 
different J-dependence from these "monopole-corrected" 
interactions GXPFM and KB3G, especially for the J = 4 
and J = 5 matrix elements. The former is more attrac- 
tive and the latter more repulsive by about 0.3 MeV. Note 
that these differences keep the monopole centroid quite 
similar value for these interactions as shown in Figs. [21 
andU 

In order to examine to what extent such differences af- 
fect the energy spectra, similar calculations have been 
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panel) 2f in ^®Ni and (lower panel) 5/2^ in ^^Ni between 
the experimental data and the results of various effective in- 
teractions. The shell-model calculations were carried out in 
the t = 7 and 6 truncated subspaces for ^^Ni and ^'^Ni, re- 
spectively. 
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FIG. 26: Comparison of two-body matrix elements 
\^(7373; Jl) between various effective interactions. 



carried for both GXPFM and KB3G interactions, by 
modifying the J = 4 matrix element to be more at- 
tractive by 0.300 MeV and J = 5 to be more repul- 
sive by 0.245 MeV so that the monopole centroid is 
kept unchanged. The modified interactions are denoted 
as GXPFM' and KB3G', respectively. The results are 
shown in Fig. It is clearly seen that the description 
is improved in both GXPFM' and KB3G', although it is 
not enough for the latter. This fact highlights the im- 
portance of modifications of the microscopic interaction 
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that go beyond the monopole type. 

The above J-dependent modifications in GXPFl 
can be understood to a good extent in terms of 
a quadrupole-quadrupole interaction [(/7/2)Up3/2)]''^'' ■ 
[(^3/2)^/7/2)]^^^- The matrix elements of this interaction 
are shown in Fig. 1261 which give similar J-dependence 
to that of the difference between GXPFl and GXPFM. 
Indeed, the y(7373; Jl)'s can be obtained by recoupling 
this interaction with the strength -0.81 MeV for GXPFl 
together with other multipoles with minor contributions, 
while this strength turns out to be —0.32, —0.34, and 
-0.38 MeV for GXPFM, KB3G and G, respectively. This 
quadrupole correction can arise from a core-polarization 
diagram where the external lines are /7/2 and with 
exchange between them and the bubble of J'" = 2+ is 
involved. If this is the case, the coupling should be par- 
ticularly strong in this channel and/or the energy denom- 
inator should smaller, as compared to what is assumed in 
the microscopic calculation of effective interaction. This 
is an intriguing problem of the effective interaction, and 
further studies are needed. 



IV. SUMMARY 

Summarizing this paper, the effective interaction for 
p/-shell nuclei, GXPFl, has been examined and tested 
from various viewpoints. It was obtained by adding em- 
pirical corrections to the microscopic interaction derived 
from the nucleon-nucleon potential. The corrections were 
determined through the systematic fitting to experimen- 
tal energy data. The most significant modifications are 
those for the monopole and pairing parts of the Hamilto- 
nian The monopole plus pairing corrections give a sub- 
stantial improvement, but are not sufficient for the de- 
scription core-excitations over the N ox Z — 28 shell gap. 
Other modifications are necessary including some which 
can be related to the quadrupole-quadrupole component. 

The analysis of the ground-state wave functions ob- 
tained by GXPFl shows that the assumption of an inert 
^^Ni core is not realistic even for the Ni isotopes. Thus, 
the amount of core-excitations strongly depends on both 
valence proton and neutron numbers relative to the core. 

The calculated binding energies agree with the exper- 
imental data quite well over a wide mass range even for 
many nuclei which were not included in the fit. How- 
ever, in neutron-rich nuclei {N > 34) around Z ^--^24, 
the difference between theory and experiment becomes 
larger, suggesting a need for including the gg/2 orbit for 
describing large deformation. Such a discrepancy can 
also be seen in the calculated 2^ excitation energy for 
^'^Cr, which is predicted too high in comparison to the 
recent experimental data. 

The description of the magnetic dipole and electric 
quadrupole moments is basically successful with a few 
exceptions. The most remarkable difference between the- 
ory and experiment can be found in the 2f states of Zn 



isotopes, which also suggests influence of the 519/2 orbit. 
As for the electric quadrupole moments, more precise 
experimental data are desired to discuss the quality of 
description especially for Ni and Zn isotopes. 

The energy spectra and electro-magnetic transition 
matrix elements for iV or Z = 28 semi-magic nuclei have 
been calculated and compared with experimental data. 
The effects of the core-excitation across N, Z = 28 shell 
gap can be seen already in the non-yrast states around 
Ex ^ 3 MeV as well as high spin yrast states. The present 
results show the reliability of GXPFl even in the cases 
that the core-excitation plays a crucial role. Consider- 
ing the current experimental and theoretical situations, 
it would be difficult to explore such significantly core- 
excited states in other nuclei with more "valence" parti- 
cles or holes relative to ^^Ni core, due to the explosive 
increase of the low-lying level density. Therefore the de- 
scription of ^^Ni and neighboring nuclei is an important 
test of the effective interaction. 

It has been shown that the reasonable description of 
N = Z odd-odd nuclei can be obtained by GXPFl, indi- 
cating its proper isospin dependent structure, although 
the experimental data are limited especially for the tran- 
sitions. 

Most of shell-model calculations were done in the 
conventional method so as to reach higher accuracy. 
The truncations were carefully examined to be accurate 
enough and for the data we compare with. In some 
shape coexistence cases, strongly deformed states appear 
at lower energies, where the Monte-Carlo Shell Model 
still produces more accurate results as compared to large- 
scale conventional calculations with high (i.e., unrestric- 
tive) truncations. Thus, one has to be cautious in using 
conventional calculations for such cases. This may be 
an important lesson for further studies on heavier nuclei 
with more developed collectivity. 

Further applications to unexplored regimes of large 
proton/neutron numbers or high excitation energy is of 
great interest. Future experiments will test the predic- 
tions and provide guidance for further improvements in 
the interaction. 
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